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Several filtering methods for the detection of gravitational wave bursts in interferometric 
detectors are presented. These are simple and fast methods which can act as online 
triggers. All methods are compared to matched filtering with the help of a figure of 
merit based on the detection of supernovae signals simulated by Zwerger and Miiller. 

1. Introduction 

Supernovae have been historicaUy the first envisaged sources of gravitational waves 
(GW). Although binary inspirals or even periodic GW emitters like pulsars seem to 
be nowadays more promising sources, impulsive sources of GW such as supernovae 
should also be considered in the data analysis design of interferometric detectors 
currently under construction (LIGO, VIRGO). Impulsive GW sources are typically 
collapses of massive stars, leading to the birth of a neutron star (type II supernova) 
0'i'i or of a black hole □; mergers of compact binaries can also be considered as 
impulsive sources i. 

The problem with such sources is that the emitted waveforms are very poorly 
predicted, unlike the binary inspirals. As a consequence, this forbids the use of 
matched filtering for the detection of GW bursts. The filtering of such bursts 
should therefore be as general and robust as possible and with minimal a priori 
assumptions on the waveforms. A drawback is of course that such filters will be 
sensitive to non-stationary noise as well as to GW bursts; spurious events, e.g. 
generated by transient noise, should be eliminated afterwards when working in 
coincidence with other detectors. But, on the other hand, burst filters could help 
to identify and understand these noise sources, which would be useful especially 
during the commissioning phase of the detector. 

All the filters presented here are dedicated to GW bursts detection and are 
compared by studying their performance to detect a reference sample of GW burst 
signals, numerically computed by Zwerger and Miiller (ZM).H 

Throughout the following, we assume that the detector noise is white, stationary 
and Gaussian with zero mean. For numerical estimates, we chose the flat (ampli- 
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tude) spectral density to be ft,„ ~ 4 x 10^^'^/\/Hz and the sampling frequency 
fs ~ 20 kHz, so the standard deviation of the noise is cr„ = hn\/ ~ 4 x 10"^"'^. 
The value chosen for /i„ correspoiids approximately to the minimum of the sensitiv- 
ity curve of the VIRGO detector around this minimum, the sensitivity is rather 
flat, in the range [200 Hz,lkHz], which is precisely the range of interest for the 
gravitational wave bursts we are interested in. This validates then our assumption 
of a white noise ; otherwise, we can always assume that the detector output has 
been first whitened by a suitable filter EEl. 



2. General filters 



2.1. Filters based on the autocorrelation 

The noise being whitened, its autocorrelation is ideally a Dirac function and in 
practice vanishes outside of zero. The autocorrelation of the data x(t) 

A^{t) ^ j x{t)x{t + T)dt (1) 

should then reveal the presence of some signal (which is surely correlated). The 
information contained in the autocorrelation function can be extracted in different 
ways. We have studied two of them and built so two non-linear filters. Jhe first 
one computes the maximum of Ax (t) and has already been described in In the 
following, we will refer to this filter as the Norm Filter (NF). A similar approach 
has been developed independently by Flanagan and Hughes in the context of the 
detection of binary black hole mergers llil. 

Another possibility is to look at the norm of the autocorrelation function (NA 
Filter) : 



ll^ll 
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]^E^W^' (2) 
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where A{k) denotes the discrete autocorrelation of N data Xi. The sum is here 
initiated at the second bin according to the fact that the noise (uncorrelated) con- 
tributes essentially to the first bin. Note that the only parameter for these two 
filters is the window size N . The behavior of the NA filter with noise only is not 
known analytically and its characteristics (mean and standard deviation) have to 
be found numerically (adding some complexity to this filtering method). 

2.2. The Bin Counting method 

This filter (BC) computes the number of bins in a window of size N whose value 
exceeds some threshold s x (t„. The threshold s is chosen by maximizing the signal 
to noise ratio (SNR) when detecting the signals of the ZM catalogue (for more 
details, see u). 



2.3. Linear Fit Filters 
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This filter fits the data to a straight line in a window of size N. If the data are 
pure white noise with zero mean the slope and the offset of the fitted line are 
zero on average, so this can well discriminate between the two cases : only noise 
or noise+signal. The slope and the offset of the fit can easily be computed as a 
function of time and of the Xi (see 0) . 

In fact, the slope and the offset are two correlated random variables. By comput- 
ing their Covariance Matrix, one reduces them to two uncorrelated normal random 
variables, X+ and with 

where X = SNR(slope) = \slope\/ asiope, Y = SNR(offset). Finally, the sum 
X"^ + X'^ is a like random variable and gives us a new filter, ALF (Advanced 
Linearfit Filter). 

2.4. The Peak Correlator 

Filtering by correlating the data with peak (or pulse) templates is justified by the 
fact that simulated supernovae GW signals exhibit one (or more) peaks. The pulse 
templates have been built from truncated Gaussian functions. The method and 
results are explained in i. 

3. Performance and efficiency of the filters 

3.1. Definition of a false alarm rate 

We arbitrarily set the false alarm rate for each of the filters to be 10~^ (72 false 
alarms per hour for a sampling frequency fs = 20 kHz). This high rate is required 
because the signals we look for are very weak. False alarms will be discarded later 
when working in coincidence. 

3.2. The Zwerger and Miiller Catalogue 

The catalogue of Zwerger and Miiller i contains 78 gravitational- wave signals. Each 
of them corresponds to a particular set of parameters (e.g initial distribution of 
angular momentum). All the signals are computed for a source located at lOMpc. 
We can then re-scale the waveforms in order to locate the source at any distance d. 

Since the signal waveforms are here known, we can explicitly derive the optimal 
SNR provided by the Wiener filter matched to each of them, and then compute the 
maximal distance of detection. We will then be able to build a benchmark for the 
different filters by comparing their results (detection distances) to the results of the 
Wiener filter (we consider here optimally polarized GW's, along the interferometer 
arms) . 
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The mean distance obtained for the Wiener Filter, averaged over all the signals, 
is about c?opt — 25.4 kpc, which is of the order of the diameter of the Milky Way. 

3.3. Estimating a filter power 

The optimal (Wiener) filtering allows to detect the i"^ signal in the Catalogue emit- 
ted by a source located up to a distance d/^^ . Similarly, a filter T is able to detect 
the same signal up to a distance c?/'^^ ; of course dj^'^^ is averaged over many noise 
realizations (about 1000) in a Monte Carlo simulation. The detection performance 
of the filter T for this signal is simply defined as the distance of detection relative 
to the optimal distance of detection : dj^'^'^ / d^"^^ . The global performance of T is 
then estimated as the detection performance averaged over all the waveforms of the 
catalogue: 

1 ^® d-^-^^ 

i=l "i 

For a given filter, and a given source located at a distance d, one can also 
evaluate a detection efficiency e, which is the number of detections n over the total 
number of noise realisations A/". This efficiency (averaged over all the signals of the 
catalogue) will characterize the practical behaviour of the filter. 

3.4. Comparison of the filtering methods : Performance 

The results for the different filters are reported in Table [l]. We also give the average 
distance of detection d= Y^'iLi '^'^^ for all the filters. 

Table 1: Performance of the different filters. L means linear filter and NL means 
non-linear filter. 
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L 


NL 


NL 


NL 


L 


L 


L 
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The three first filters NF, NA and BC (all non-linear) have a performance slightly 
below one half, while the PC have a performance greater than 0.7. Both X+ (or 
Ar_) and ALF can reach a performance around 0.8. Note that ALF has been in fact 
implemented with a sampling of 20 different window sizes, sufficient to cover the 
variety of signals (with a non-significant loss of generality). If implemented with a 
single window size, as the other filters NF, NA and BC, its performance decreases 
down to 0.72. 

3.5. Comparison of the filtering methods : Detection efficiency 

For the Wiener filter, one can show that the mean efficiency is roughly 50% for 
signals located at a distance dl'^\ If p is the mean performance of a given filter. 
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one would expect a 50% detection efficiency at ~ p x . In fact, such 

efficiency is reached for a smaller distance. 

If one defines the effective performance p^f / as the ratio dj^"^'' /d^'^^ for which the 
detection efficiency is about 50%, then peff — 0.74 for ALF and X_ and p^ff ~ 0.71 
for X^. This definition gives an idea of the efficiency one can reach in practice, and 
has to be taken into account when choosing between different online triggers. 

4. Conclusion 

We have discussed several filters to be used as triggers for detecting GW bursts in 
interferometric detectors. They are all sub-optimal but their performance is close 
to the one obtained with the Wiener Filter. 

Concerning the detection of Zwerger-Miiller-like signals, we note that none of 
the BC, NF and NA filters is efficient enough to cover the whole Galaxy in average 
(but their window sizes have not yet been optimized), contrary to ALF and PC 
(and optimal) filters. A few signals can be detected at distances beyond 50kpc, the 
distance of the Large Magellanic Cloud (LMC). It is clear that this class of signals 
will be detected by the first generation interferometric detectors such as VIRGO 
only if the supernovae occur inside our Galaxy or in the very close neighbourhood. 

Finally, all the filters studied here can be implemented on line without problem, 
due to use of FFT's (for the NA and the PC) or to simple recursive relations between 
filter outputs in successive windows (NF,BC or ALF). 

Correlations and coincidences between those filters are under study in order to 
either reduce background (hence a quantifiable loss of signal) or lower detection 
thresholds (hence a gain of a few % in performance and efficiency). 
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